Macroscopic Periodic Tunneling of Fermi Atoms in the BCS-BEC Crossover 
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We study the macroscopic quantum tunneling of two weakly-linked superfluids made of interacting 
fermionic atoms. We derive atomic Josephson junction equations and find that zero-mode and n- 
mode frequencies of coherent atomic oscillations depend on the tunneling coefficient and the sound 
velocity of the superfluid. By considering a superfluid of 40 K atoms, we calculate these oscillation 
frequencies in the crossover from the Bardeen-Cooper-Schrieffer state of weakly-bound Cooper pairs 
to the Bose-Einstein Condensate of strongly-bound molecular dimers. 
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I. INTRODUCTION 

The crossover from the Bardeen-Cooper-Schrieffer 
(BCS) state of Cooper fermion pairs to the Bose-Einstein 
condensate (BEC) of molecular dimers with ultra-cold 
two-hyperfine-cornponents Fermi vapors of 40 K atoms 
[ESS] and 6 Li atoms [3, H, has been observed in 
the last years by several experimental groups with the 
use of Fano-Feshbach resonances Q ■ From the theoreti- 
cal side, extended BCS (EBCS) equations [1,1^, lioj have 
been used to reproduce density profiles and collec- 
tive oscillations [T3. [l3l| of these Fermi gases. In addition, 
more recent calculations based on Monte Carlo (MC) fit- 
ting and superfluid dynamics [l4| have shown that the 
mean-field EBCS theory is quite accurate. 

Recently, Spuntarelli et al. [HJ studied the stationary 
Josephson effect [l6[ across the BCS-BEC crossover with 
neutral fermions by using the EBCS equations: in detail, 
they computed the current-phase relation throughout the 
BCS-BEC crossover at zero temperature for a two-spin 
component Fermi gas in the presence of a barrier. The 
Josephson effect of atomic superfluids, i.e. coherent os- 
cillations between two weakly- linked bosonic clouds, was 
predicted [ItJ and observed [l8| with BECs. Josephson 
oscillations in superfluid atomic Fermi gases have been 
thoretically considered by Paraoanu et al. [Hj], Wouters 
et al. 20], and Adhikari [2lJ. The macroscopic oscilla- 
tions of tunneling neutral atoms is closely related to the 
familiar Josephson effect of charged electrons in super- 
conductor junctions [22j and its investigation can help a 
deeper understanding of these exotic states of matter. 

In this work we present an investigation of the dynam- 
ical Josephson effect in the BCS-BEC crossover based 
on a time dependent local density approximation for the 
dynamics of the Ginzburg-Landau (GL) order parameter 
of the atomic Cooper pairs [23| , starting from a reliable 
parameterization of the bulk chemical potential [l4| in 
the crossover. From our zero-temperature GL equation 
we obtain the atomic Josephson junction equations for 
two weekly linked fermionic superfluids. The main result 
provides zero-mode and 7r-mode Josephson frequencies 



of periodic quantum tunneling as a function of the sound 
velocity. Our investigation of the Josephson effect with 
neutral Fermi atoms, a direct manifestation of macro- 
scopic quantum phase coherence, is of interest not only 
conceptually but also for future applications in quantum 
computing [H HI]. 



II. GINZBURG-LANDAU EQUATION AT ZERO 
TEMPERATURE 



In a dilute Fermi gas of N atoms with two equally 
populated spin components and attractive inter-atomic 
strength at zero temperature, superfluidity and coher- 
ence are strongly related to the properties of the two- 
particle density matrix. Let Vv(r, t) be the field operator 
that destroys a fermion of spin a (a =T,i) in the posi- 
tion r at time t. The two-particle density matrix can be 
written as 



(n, t)$+ (r 2 , t)Vv ; (r' 2 ,t)yl a , (r' 2 , t)) = 
N J ( riCri ' r2<72 > QXi ( r i CT 'i ' r 2°2 . *) 



(1) 



where ( • ) is the ground-state average, Xj an d Nj are 
respectively normalized eigenfunctions and eigenvalues 
of the two-particle density matrix, and the maximum 
eigenvalue cannot exceed N/2 [3. [2q|. If one eigenvalue, 
say Nq, is of order N/2, off-diagonal long-range order is 
present. 

At zero temperature, where the superfluid density co- 
incides with the total density, the GL order parameter 
describing the motion of Cooper pairs of atoms is defined 
as 



*(r,t) 



Xo (rT,r|,t) = y^Mex P (^(r,t)), 

(2) 

where n(r, t) is the local atomic number density (n(r, t) /2 
is the local density of pairs) and 8(r,t) is precisely the 
phase of the condensate wave function xo l23l H3]- 
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Note that 9(r,i) is also the phase of the gap func- 
tion A(r, t) = |A(r, t)\ exp (i9(r,t)) of Cooper pairs 
0, US EE E3] ■ Under an external potential U (r) acting 
on individual atoms, the low-energy collective properties 
of the Fermi superfluid [lj] can be described by the the 
following highly nonlinear time-dependent GL equation 
(TDGLE) 



injL*(r,i) 



4m 



- 2 J7(r) + 2 /*(n(r)) 



*(r,t). 



Here m is the mass of one atom and 

d 



(3) 



(4) 



is the atomic bulk chemical potential of a homogeneous 
fluid with density n 28], if £(n) is its energy per particle. 
The phase of the order parameter drives the superfluid 
velocity 



v(r,*) = _V*(r,t), 



(5) 



which is irrotational (V A v — 0) by construction [28j |. 
The equation of superfluid velocity permits to map Eq. 
(|3|) to the hydrodynamic equations of fermionic superflu- 
ids: 



On 



V • (nv) = , 



<9v „ 
m— + V 

ot 



H 2 W 2 Jn m 2 , ' 



(6) 
0.(7) 



In these zero-temperature hydrodynamical equations, 
which are valid in the full BCS-BEC crossover to de- 
scribe macroscopic long- wavelength phenomena of the 
fermionic superfluid [29j . statistics enters through the 
equation of state fJ,(n) and the quantum-pressure term 
— [h 2 / (8my/n)]V 2 y/n, which is absent in the classical hy- 
drodynamic equations [H, [2i|. The coefficient 1/8 in 
the quantum pressure holds for pairs of fermions, as op- 
posed to atomic bosons, where the coefficient is 1/2, and 
v = (h/m)W9. The bulk chemical potential fi(n) is the 
key ingredient of our study; theoretical calculations in 
both asymptotic limits of 1/y and Montecarlo data (30j 
suggest that in the BCS-BEC crossover it can be written 



with ep — (ti 2 k 2 7 /2m) the Fermi energy and kp — 
(37r 2 n) 1/3 the Fermi wave number. The parametrization 
of e(y) chosen in Ref. [l4| is the following 



e(y) = ai — «2 arctan 013 y 



Pi 



Ih 



1/7 1 



(10) 



where the values of the parameters ai, a-2, a$, j3i, 0% are 
reported in that work. In the BEC regime (i/> 1), where 
fx(n) ~ n, the TDGLE reduces to the time-dependent 
Gross-Pitaevskii equation for composite bosons of mass 
2m subject to the effective potential 2E/(r) + 2/i(n(r, t)) 
[TH - In the BCS regime (y -C —1), where fj,(n) ~ 
€f {1 + kpdF / (3tt)), the TDGLE gives the Anderson- 
Bogoliubov mode (sound velocity) of neutral supercon- 
ductors. 

Some years ago an effective nonlinear Schrodinger 
equation (NLSE) for superconductors was derived from 
the microscopic BCS Lagrangian in the low-frequency 
long-wavelength limit of |A(r,t)| by Aitchison et al. and 
De Palo et al. 131]. That NLSE and our TDGLE pro- 
duce the same Anderson-Bogoliubov mode [3, [3l| , but 
TDGLE describes accurately sound velocity and other 
collective modes of the fermionic superfluid in the full 
BCS-BEC crossover [13. 



III. DYNAMICS OF TWO WEAKLY-LINKED 
FERMI SUPERFLUIDS 

Assume that ?7(r) is a potential with a barrier that 
splits the superfluid into two subsystems A and B sep- 
arated by a region C where the modulus of the order 
parameter is exponentially small [32] . We look for a time- 
dependent solution of the TDGLE of the form 



tf (r,t) = * A (t) * A (r) + *bW *n(r) , 



(11) 



where <fr a (r) is the quasi-stationary solution (real and 
normalized to unity) of the TDGLE localized in region a 
(a = A, B). 

By inserting this ansatz for \? into Eq. ([3]), after in- 
tegration over space and neglecting exponentially small 
terms, the system can be described by the follow- 
ing two-state model (33| : 



,W^(3, 2 n) 2 / 3 



(8) 



where e(y) is a dimensionless universal function of the 
inverse interaction parameter y — 1/kpap, where kp — 
(37r 2 n) -1 / 3 is the Fermi wavenumber of the noninteract- 
ing fermions and of is the fermion-fermion scattering 
length. The function e(y) was parameterized across the 
BCS-BEC crossover by Manini and Salasnich [l4| to fit 
the MC data of Astrakharchik et al. [30] and the asymp- 
totic expressions of the bulk energy per particle 



r£Fe(y) 



(9) 



ot 
d 



E A * 



A *A 



K * T 



at 

for the two complex coefficients ^ a (t). Here E a 
is the energy in region a, with 



E 



4m 



2U 



and 



Ei ~ y"$ a 2/i(2|* Q | 2 $ 2 )$ Q d 3 



(12) 
(13) 

K + 

(14) 
(15) 
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The coupling energy K describes phenomenologically the 
tunneling between the two regions. It is quite common to 
use a phenomenological tunneling energy to analyze the 
Josephson effect in superconductors M, I22I I27I [H| . In 
general, thermal [3~i| and quantum [351 ] fluctuations will 
affect the value of K. 

At zero-temperature the single-particle tunneling en- 
ergy to can be estimated as 



to 



4m 



2U 



<f> B d 3 r, (16) 



if the barrier is orthogonal to the x axis and such that 

fc iw ~ Kn))] 1/2 dx > 1. In the BEC region 

the quantum depletion is negligible [26j . the tunneling is 
fully coherent, and the coupling energy is simply K ~ 2to 
[13, [HI]. Instead in the deep BCS regime, where there is 
a very large quantum depletion, only a small fraction of 
pairs perform coherent tunneling, and microscopic cal- 
culations suggest K ~ |A|ig/(87re|7V) [|f|. Unfortu- 
nately a microscopic derivation of if in the full BCS-BEC 
crossover is not yet available. 

Under the assumption that the double-well potential 
U(r) has the shape of two weakly- linked sharp-edged 
boxes of volumes Va and V B , we can write *f/ a (t) = 
y/N a (t)/2exp (i0 a (t)), where N a (t) and 9 a (t) are the 
number of fermions and the phase in region a. In terms 
of the phase difference 



<p(t) = 9 B (t) - e A (t) 

and relative number imbalance 

N A (t)-N B (t) 



z(t) 



N 



(17) 



(18) 



with N = N A (t) + N B (t) the constant total number of 
atoms, Eqs. (fl~2|) and (p~3|) give 



<P = T 



2K 

2 " 
h 

2K 



( N , 



cos ip 



E° A -E B 



(19) 



(20) 



These are the atomic Josephson junction (AJJ) equations 
for the two dynamical variables z(t) and <p(t) describ- 
ing the oscillations of N Fermi atoms tunneling in the 
superfluid state between region A of volume V A and re- 
gion B of volume Vb. These equations, linking the tun- 
neling current I = —zN/2 = (KN/H) yl — z 2 sia(p = 
Io y \ — z 2 sin tp to the phase difference tp, reduce to the 
familiar Josephson's expression I = Iq sin(^), well estab- 
lished for BCS superconductors, in the appropriate limit 
\z\ <C 1 dj|. Also, in the deep BEC regime, where fi(n) ~ 
n, the AJJ equations reduce to the bosonic Josephson 
junction (BJJ) equations introduced by Smerzi et al. (l7| 
and found to describe accurately the density oscillations 
of weakly linked condensates [18| • 



IV. AJJ EQUATIONS: LINEAR REGIME 

The AJJ equations can be solved across the BCS-BEC 
crossover, even taking into account a small but finite im- 
balance z. To first order in z, Eqs. (fill)) and (|2U|) read 



2K 

— ~h sm v ' 

2A 2K 

— z H — —z cos tp - 

h h 



2Aab 



where 



and 



A= — 



N 


' 1 dfj, 




1 9/x 




~2 




H 

A 




B 



MAB 



-( M | A - M | B ) + (E° A -El) 



(21) 
(22) 

(23) 

(24) 



with 



Ml 



and 2«i 

on 



the bulk chemical potential and its 

derivative are calculated at the density N/(2V a ), a = 
A,B. 

In the symmetric case (V A = V B = V/2 and E° A = E^) 
A takes the particularly simple form 



A = 2n 



dn 



2mc a 



(25) 



where c s is the sound velocity computed at the mean 
density n = N/V of the superfluid. ^,From Eq. (JSJ one 
has thus 



A = 2mvp 



i e ( y ) _ l e '( y ) + y- e "(y) 
3 yy> 5 Ky> 30 Ky ' 



(26) 



where vp = hkp/m is the Fermi velocity of non- 
interacting fermions. 

In the symmetric case jiAB = 0, Eqs. (j2Tj) and (|22|) 
admit a stable stationary solution (z, Cp) with z — and 
<p = 2irj, for integer j. If A/K < 1, also z — and 
93 = 7r(2j + 1) is a stable stationary solution. These 
stationary solutions remain stable even if one includes 
higher-order terms in the z-expansion of Eqs. (|19p and 
(f2"0")) . Small oscillations around the stable solutions with 
Cp — 2irj and (p = n(2j + 1) have frequencies 



K 



A 



(27) 



and are called zero-mode (with +) and 7r-mode (with — ), 
respectively. The zero-mode is the analog of the Joseph- 
son plasma frequency in superconducting junctions [22|. 
The analogous of this 7r-mode was observed in weakly 
coupled reservoirs of 3 He-B [37| and discussed in the BJJ 
equations [l7| . 

By using the functional dependence on the inverse in- 
teraction parameter y = \/(kpa,p) of the sound velocity 
c s , and inserting this expression of A in Eq. (|27[) . we ob- 
tain the oscillation frequencies vo and v^. As an example, 
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- iJ=l|im y=i.i9 v=4i.3 Hz 

,'\ >\ A ,'\ /\ r< 



1 2 

y = l/(k F a F ) 



FIG. 1: (color online). Zero-mode frequency vo (solid) and n- 
mode frequency v v (dashed) for a superfluid of N = 10 6 40 K 
atoms between two symmetric regions of volume 25 • 10 6 /im 3 
(mean density n = 0.02 /im -3 ) and fixed tunneling parameter 



K/ki 



10" a Kelvin. 



consider a fermionic superfluid of K atoms with total 
density n — 0.02 atoms//im 3 : Fig. [T] reports v / n as a 
function of y, assuming that the tunneling energy can 
be held fixed to the value K/kn — 10~ 9 Kelvin (where 
fee is the Boltzmann constant). As previously stressed 
in the deep BCS regime K is proportional to |A| and 
becomes rapidly very small, thus in practice it would be 
difficult to keep constant. The zero- mode frequency de- 
creases for increasing y across the crossover, approaching 
the asymptotic value K/(ttH) (dotted line) in the BEC 
limit (i/ » 1). In contrast, the 7r-mode does not exist 
until past the unitarity point y = 0, i.e. until the value of 
A becomes smaller than K; then increases approach- 
ing the same large- y limit as isq. 

In limiting cases of the BCS-BEC crossover analytical 
expressions are available for A. (i) Deep BCS regime 
(y <C —1): the sound velocity c s approaches vf/VS 
PJ, and A = 2mv% = 2(/i 2 /to)(3tt 2 n) 2 / 3 . (ii) Unitar- 
ity point (y = 0): according to MC simulations of As- 
trakharchik et al. (30j the sound velocity is c s ~ 0.37 Vp, 
thus A = 0.29 mv 2 ^. (iii) Deep BEC regime: all molec- 
ular dimers are in the BEC and the sound velocity is 
such that mc 2 s — nh 2 a]\fn/m where au = 0.6 Of 
is the dimer-dimer scattering length 30]. In this limit 
ap — * + for non- interacting bosonic dimers, A = and 
corresponding vq = V v = K/ (nh) . 



V. AJJ EQUATIONS: NONLINEAR EFFECTS 

We come now to nonlinear effects in the AJJ equations 
(fig)) and d20D in tue symmetric case (V A = V B = V/2 
). The variables if and z are canonically 




FIG. 2: (color online). Nonlinear effects on the zero-mode 
oscillation of N = 10 6 40 K atoms between two symmetric 
regions with the same conditions of Fig. [T] at fixed interaction 
clf = 1 /im, corresponding to y — 1.19. These trajectories 
result from the integration of the coupled equations (|19p and 
(|20[) for z(t) (solid) and ip(t) (dashed), with initial conditions 
<p(0) = 0, and z(0) = 0.5 (left) or z(0) = 0.999 (right). 



conjugate, with 



9 



where the Hamiltonian is 



dH 

dip 
dH 
dz 



2K 



z* cos ip 



(28) 
(29) 

(30) 



and E\ 



E B) 



with 

G'(z) = F(z) = fi(y (1 + z)) -J?L(1- Z )j . (31) 

It is straightforward to show that these nonlinear equa- 
tions admit a symmetry-breaking stationary solution 
(z,(p), with Cp — Tr(2j + 1) and with z the solution of 
the equation F(z) = K zj\J\ — z 2 . The analysis of sta- 
bility shows that this symmetry-breaking solution, where 
the superfluid displays macroscopic self-trapping, is sta- 
ble only if K > F'(z)(l - z 2 fl 2 . 

To analyze nonlinear effects we solve numerically the 
AJJ equations. Regular zero- mode oscillations of z(t) 
and ip(t) are displayed in Fig. ® under the same condi- 
tions as for Fig. [TJ The oscillation starting from z(0) = 
0.5 indicates that the solution (|27|) of the linearized equa- 
tions (|2"Tj) and (|2"2"|) are fairly accurate even for finite and 
not quite small amplitude. Eventually however, for very 
large amplitude, z(0) = 0.999, deviations from the har- 
monic approximation become quite visible. These devia- 
tions are illustrated, for increasing amplitude, in Fig. [3] 
the left panel shows that the frequencies of the zero- and 
7r- modes approach each other. The right panel shows 
the decay in the fractional leading Fourier component of 
the oscillation <p(t), as the wave shape distorts from a 
perfect sinusoid, and acquires higher (odd) Fourier com- 
ponents. Note the extremely small deformations from 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 
z(0) z(0) 



FIG. 3: (color online). Nonlinear effects on the zero-mode 
(solid) and 7r-mode (dashed) frequency (left) and waveform 
(right) for the same conditions as in Fig. [2] as functions of 
the initial amplitude z(0). In an expansion ip(t) — ip(0) + 
S^Li a k sin(27Ti/fci), the fractional leading Fourier component 



of the right panel is a 1 / a k 



BEC crossover by using AJJ equations. We remark that 
the coupling energy K appearing in our Josephson equa- 
tions is a phenomenological parameter: from the experi- 
mental measurement of the frequencies of periodic quan- 
tum tuneling one can infer the value of K by using Eq. 
(|2"T)) . Analytical expression of K based on microscopic 
theory are available only in the deep BCS regime and in 
the deep BEC regime. An important issue is surely the 
development of a microscopic theory of tunneling in the 
full BEC-BEC crossover. We also stress that all Joseph- 
son oscillatory frequencies discussed here cannot exceed 
the frequency |A|/(7r/i), associated to the breaking of 
Cooper pairs. The AJJ equations can be extended to in- 
vestigate Josephson junction arrays for neutral fermionic 
atoms in optical lattices, potentially of extreme relevance 
for quantum information and quantum computing with 
ultracold atoms. 



perfect harmonic oscillations and the minor deviations 
of the frequencies from the linear values vq = 53.0 Hz 
and u n = 25.8 Hz of Eq. ([27]), until z(0) ~ 0.6. Sim- 
ilar evolutions near the y — point on the BEC side, 
and on the BCS side, show oscillations of z(t) around 
the self-trapping symmetry-breaking stationary solution 
z, accompanied by monotonously increasing phase f(t). 

VI. CONCLUSIONS 

We have investigated the macroscopic quantum tunnel- 
ing of two weakly-linked Fermi superfiuids in the BCS- 
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